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Abstract 

Static solutions of white dwarfs with spherical symmetry and local anisotropy 
are studied in the post-Newtonian approximation. It is argued that the condition 
for equilibrium must be that the total energy is a minimum for given baryon 
number and the question whether there is local isotropy or anisotropy inside the 
star should follow from that condition, rather than be postulated "a priori". It is 
shown show that, in post-Newtonian gravity, there are stable configurations 
with local anisotropy for masses above the Chandrasekhar limit. 
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1. INTRODUCTION 



It is well known that both Newtonian gravity and general relativity allow 
spherically symmetric solutions with local anisotropy and that anisotropy may 
have an important effect on the stability [1]. Einstein himself considered in 
1939 a model with local anisotropy, consisting of a spherical cluster of 
particles each moving in a circle, and concluded that "the essential result of this 
investigation is a clear understanding as to why the 'Scwartzchild singularities' 
do not exist in physical reality" [2]. (It is ironic that the same year 
Oppenheimer and Snyder [3] published the first calculation demostrating the 
formation of a black hole, and Einstein's paper was cast into oblivion). 

In more recent times, local isotropy has been the common assumption 
in studies of stellar evolution. Nevertheless, different scenarios have been 
proposed which incorporate local anisotropy into the modelling of stellar 
objects, e.g. pion condensation, phase transitions, boson stars and crystalization 
in white dwarfs, among others (see [1] and references therein). A standard 
assumption in these approaches is that the departure from local isotropy is 
caused by the equation of state of the matter forming the star being essentially 
anisotropic. In this paper we shall make a different assumption, namely that 
anisotropy may be produced by a dynamical effect of the gravitational field 
itself even if the equation of state of the fluid forming the star is isotropic in 
free space. 

We shall argue that the condition of equilibrium of white dwarfs is that 
the total energy of the star is a minimum for given number of baryons. We 
shall show that local anisotropy is a straightforward consequence of this 
general condition in the case of massive white dwarfs. We consider only 
Newtonian and post-Newtonian approximations, leaving for the future the 
study of stars where the gravitational field must be treated using general 
relativity. 



2. NEWTONIAN WHITE DWARFS WITH LOCAL ANISOTROPY 



A white dwarf consists of a plasma of protons, heavier nuclei (mainly 
helium) and electrons. We shall assume that the chemical composition of the 
plasma is homogeneous throughout the star. The temperature of the plasma 
inside the star may be taken as zero Kelvin (this is because the quantum 
zeropoint energy of the electrons is much greater than the thermal energy of the 
plasma, see e.g. Ref.4). As a consequence the star is, to a good approximation, 
an isolated physical system at zero temperature and the standard condition for 
stable equilibrium holds, namely the total energy should be a minimum for 
given baryon number (or, in Newtonian gravity, minimum energy for a given 
mass, because mass is proportional to baryon number in this case). We should 
not impose additional constraints. In particular we must not assume "a priori" 
that the pressure is isotropic at every point inside the star, but the general 
condition of equilibrium (minimal energy) is what will determine whether the 
pressure is isotropic or not. 

In Newtonian gravity the statement of the problem is straightforward. 
We must find the ground state of a quantum system with Hamiltonian 

H = Ej V(mj^ c^ + pj^ c^ ) + Lij (ZiZj e^ - G m; mj ) x-,-^ , (2. 1) 

where i, j label two arbitrary constituent particles (nuclei and/or electrons). Zje 
(mj , pj) is the charge (mass, momentum) of particle j , and G (c) is Newton's 
constant (the velocity of light). In the Hamiltonian (2.1) we have ignored 
magnetic interactions, electron-positron pair creation and the possibility of 
neutronization (the reaction of an electron and a proton to give a neutron and a 
neutrino). To a good approximation magnetic interactions and pair creation 
may be neglected. In contrast, neutronization may be relevant in some cases [4] 
but we shall ignore it for the sake of simplicity ( we shall comment on it in 
section 7). 

The ground state of the Hamiltonian (2.1) cannot be found without 
making some approximations. An obvious one derives from the weakness and 
long range character of the gravitational field (the electrostatic interaction is 
effectively short ranged because atractions and repulsions cancel each other at 
a long distance). Hence it follows that the t5^ical distance where the 
gravitational field changes inside the star is many orders of magnitude greater 



than the de Broglie wavelength of the constituent particles. This allows us to 
rewrite the Hamiltonian (2.1) approximating the last term in the form 



H = Lj V(mj^ + pj^ c^ ) + Lij ZiZj rij"' + Lj mj 4)( rj ) , 



<J>(rj) = -G<»I'|Limirij 



-1 



(2.2) 



where <!>( rj ) is the gravitational potential at rj and ^ the wavefunction of the 
plasma. Our problem is now to find the ground state of a plasma in a given 
external potential 4>, but constrained by the consistency requirement that $ is 
the potential created by the plasma itself. So stated the problem is still 
formidable because it involves the full machinery of quantum many body 
theory, but there is a good enough approximation which is described in the 
following. 

We may consider that the particles are independent, except for some 
averaged electrostatic interaction. As the velocity of nuclei is much smaller 
than that of the electrons, we may treat them as being at rest. Furthermore, we 
may replace the nuclei by a jellium of positive charge so that the net charge is 
zero, which is a standard approximation for a plasma at zero temperature. After 
that the plasma becomes a degenerate electron gas immersed in a background 
of positive charge. Such a gas is fully characterized by the size and form of the 
Fermi surface at every point of space. In fact, from it we may obtain all 
functions of position needed to solve the problem of minimizing the energy 
insuring gravitational equilibrium (i.e. the fulfillement of the second eq.(2.2) ). 
These functions are: the electron density, n, the energy density, u, and the 
stress tensor, py , which we get as follows (we shall use Planck units G = c = h 
= 1 from now on): 



n = (47c') 



Id^k, 



nu = (47f 



1 [V(me' + k') - me ] d'k , 



Pij = (471^)"' \ ki kj [A/(me' + k') ] A\ , 



(2.3) 



where the integrals should be extended to the interior of the Fermi surface. As 
a consequence of the spherical symmetry, which we will assume for the 
(nonrotating) star, there is one privileged direction at every point inside the 
star, the radial direction. Thus the Fermi surface at every point will possess 
rotational symmetry around the radial direction. This implies that there are only 
two principal stresses, namely the radial pressure, pr , and the transverse 
pressure, pt , for which the last eqs.(2.3) become eqs.(3.9) and (3.10) (see 
below). 

Now the problem of finding the ground state of the Hamiltonian (2.1) 
has been reduced to the problem of finding the size and form of the Fermi 
surface at every point inside the star in such a way that: 1) gravitational 
(hysrostatic) equilibrium exists at any point, and 2) the total energy of the star 
is a minimum for given total baryon number. In Newtonian gravity the 
hydrostatic condition for equilibrium is (we do not assume isotropic pressure 
although we do not exclude this possibility) 

dpr/dr + 2(pr-pt)/r = -mpo/r^ m(r) = 4 7t lo' po dr, (2.4) 

where we have introduced the mass density, po , defined by 



me being the mean baryon rest mass and |l the mean number of baryons per 
electron. (For the derivation of eq.(2.4) see eq.(5.3) below). The mass and 
energy of the star are given by 



In the present paper we shall consider only the case of small anisotropy, 
that is small deformation of the Fermi surface with respect to a sphere (the 
calculation in the general case is more involved and will be made elsewhere). 
Thus the Fermi surface may be characterized by just two functions of the 
radius, and we may write the mass density, po, the energy density, u, the radial 



po(r) = n mB n(r) . 



(2.5) 



M= Ipod\, 



E = |{ u-m/r } pod\. 



(2.6) 



pressure, pr. and the transverse pressure, pt , in terms of those two functions 
using eqs.(2.2) and (2.3) as is explained in the next section. After that we shall 
be able to find the functions Po(r), u(r), pr(r) and pt(r), fulfilling the hydrostatic 
equilibrium eq.(2.4) and making E a minimum for given M, eqs.(2.6). The 
practical procedure to do that is described in section 4. 

It is interesting to point out that eq.(2.4) would become the 
conventional equation of hydrostatic equilibrium if we were assuming isotropic 
pressure, that is 

dp / dr = - m po / r^ , p = pt = p^ . (2.7) 

hi this case both requirements, the energy E ( eq.(2.6)) being a minimum and 
the hydrostatic equilibrium (eq.(2.7) ) are equivalent conditions. However, if 
we do not assume local isotropy "a priori", the hydrostatic equilibrium ( now 
given by eq.(2.4)) and the minimal energy are independent conditions and both 
should hold. 

3. ANISOTROPIC DEGENERATE ELECTRON GAS. 

As said above, we shall take the electrons as noninteracting particles, although 
corrections to this simple model are well known. In a degenerate electron gas, 
anisotropy appears if the Fermi sphere of the electrons is deformed. If the 
deformation is small we may write 

kF(e) = kFo[l-yP2(z)], z = cose, (3.1) 

where is the angle with respect to a given direction (the radial direction in the 
case of a star), P2(z) is a Legendre polynomial and y is a parameter measuring 
the fractional anisotropy. In (3.1) a dipole term is excluded in order that there is 
no net electric current at any point, but higher multipoles of even order should 
be introduced if the anisotropy were not small. 

From (3.1) it is easy to get the number density of electrons (compare 
with eq.(2.3) ) 



n = (InY dz lo"^ dk = (37U^)"^ kpo^ (1+3/5 y^+2/35 y^). 

(3.2) 

The energy per electron, u, or the energy density, nu, may be also obtained 
easily. We have 

n u = (2nY 1.1^ dz L'^k^ dk [V(me^ + k^) - me ] = 

= (27l^)-^me^ dz IJ'^'^^' [V(l+ v^s^) - 1 ] ds, 

(3.3) 

where me is the electron mass and we have introduced the new variables 

V = kpo / me , s = k / kpo . 

(3.4) 

The calculation of (3.3) is straightforward using the expansion 

lo^"' f(s) ds = lo^ f(s) ds - 8 f(l) + 1/2 8^ (df/ds)s=i + 0(8^). 

(3.5) 
We get 

n u = 7U- VeS' ( X + ^ y' ) + O(y') , u = 3 me [ % + (^ - 3/5 %) y' ] + O(y'), 
(3.6) 

where we have introduced the functions 

X(v)=Jo^ [V(l+ vV) -l]ds = (SyY{y (1+ 2v^) a/(1+v^) - log[v+V(l+v^)]}, 
(3.7) 



^(v) = 1/10 (d/ds[s^ (V(l+v^s^) -1)] )s=i = 1/10 (2+3v^) /V(l+v^) - 1/5. 
(3.8) 

Similarly, we may obtain the radial and transverse pressures 
Pr = (InY J.i^ dz Jo'^k^ dk [A/(me^ + k^)]"^ = 

= (27UV)- VeV IV dz L'-''''^ s' [V(l+ ds = 

= (37c')-' me^ ( X' - 2/5 [v/V(l+ v')] y + 11/7 y' ) + O(y^) , 

(3.9) 

Pt = (2nY I.i^ 1/2(1- z^) dz lo'^k^dk [V(me^ + k^)]"^ = 

= i3nY me' ( X' + 1/5 [v/V(l+ v')] y + 5/7 y' ) + O(y^) , 

(3.10) 

where %'= d%/dv and (j)' = d(j)/dv. The mean pressure is 

p = 1/3 (pr + 2 pt ) = QtiY m,Wr + <^'y^) + O(y') . (3.11) 

It is not difficult to check that p is related to the electron density, n, and the 
energy per electron, u, by a standard thermod5niamic relation which, in the case 
of anisotropy, should be written 

p = n^ 3u/3n|y=const = (5u/3v)y /(3n/3v)y . (3.12) 

It is convenient to replace the anisotropy parameter y by another one x: 

X = (Pt - pO/3 p = 1/5 [v/V(l+ v')] y - 2/7 ^' + O(y') . (3.13) 



From the definition it follows that x € [-1, 1/2]. If we eliminate v, kpo and y 
amongst eqs.(3.2), (3.6) and (3.13) we may obtain the energy per electron, u, in 
terms of the electron density, n, and the anisotropy parameter, x: 

u = 3 me { X(v) + [25 ^(v) - 15 x(v) ] (1+v') v"' xiyY ' x^} + 0{x^) , (3. 14) 

2 1/3 

where v stands here for (37U n) /me . In particular, in the nonrelativistic 
(v«l) and ultrarelativistic (v»l) regimes we have 

u = 3/10 On\f'^lm^ (1+x^) if v«l, u = 3/4 (37C^n)^^^ (1+ 5/8 x^) if v»l. 
(3.15) 

Hence, the mean pressure, p, may be obtained from (3.12) and, taking (3.13) 
into account, the radial and transverse pressures are given by 

Pr=(l-2x)p, p, = (l+x)p. (3.16) 

Also we shall redefine u as the energy per unit mass dividing the 
previous u by |i.mB. We are interested in white dwarfs of high mass (close to 
the Chandrasekhar limit) where the electrons are ultrarelativistic. But the 
approximation of the second eq.(3.15) is not enough and we require corrections 
to that expression up to 0(v" ). Nevertheless we shall consider only cases 
where the local anisotropy is small (x«l) and this allows us to introduce 
corrections for anisotropy only in the leading term of the equation of state. 
Consequently we shall write 

u = uo+ ui+ ua , uo= 3 K po^^^ ui = K' + 3 K"po '^^ , ua = 15/8 K po^^^ x^ , 
(3.17) 

where the constants K, K'and K" may be obtained from (3.14) and (2.5). Their 
values are well known (see e.g. Ref. 4) and only the numerical coefficient of 
the anisotropy correction ua is a new result. In our work of the following 
sections both ui and ua will be considered corrections of the same order to the 



main term uq. The mean pressure may be obtained from (3.17) and the 
thermodjmamic relation (3.12). We get 

p = P0+ pi+ Pa , Po= K po^^^ pi = - K"po'^^ , Pa = 15/8 K po''' . 
(3.18) 

In a homogeneous gas without external forces, the positivity of the 
anisotropy correction (term with x in eq.(3.17)) implies that the minimal 
energy density for given density of particles corresponds to x=0, that is local 
isotropy. But we claim that the condition of minimal energy is the only cause 
why the proton-electron plasma at zero temperature is isotropic in free space. 
As is very well known the Fermi surface of the electron gas may be deformed 
by electrostatic fields (e.g. in metals) and the gas pressure is no longer 
isotropic. We think that there is no reason to postulate that it cannot be 
deformed also by the gravitational field, a postulate nevertheless made in 
conventional teatments of relativistic stars. 

4. NEWTONIAN EQUILIBRIUM WITH LOCAL ANISOTROPY 

After the results of the two previous sections, the theory of Newtonian 
white dwarfs is straightforward. (If the mass is low enough the star behaves 
like a polytrope with y = 5/3, the equilibrium state corresponding to isotropic 
pressure [5], but here we shall be concerned with the more interesting case of 
high mass stars, whose equation of state is given by (3.18).) In the conventional 
treatment [4] it is postulated that Pr = Pt , so that the problem reduces to solving 
eq.(2.7) (see the last paragraph of section 2). In practice a variational method is 
used in which a one-parameter family of configurations is selected with the star 
approximated by a polytrope. That is, we change variables: 

Po(r) = Pc e^ r = [ 47U ie'(^i)l Y'" M'^' p,-'" ^ , 

(4.1) 



and assume that 9(^) is the solution of the Lane- Emden equation of index n=3. 
The central density pc is the parameter of the family of configurations. Thus 
eq.(2.6) gives the energy of the star in terms of standard integrals and we 
obtain 

Ei = ( A M - B M^^^) pc - K' M + C M pc (4.2) 

where A, B, C and K' are well known positive constants [4] and the subindex 
of Ei stands for "isotropic". The condition of stationarity of the energy for 
given mass (i.e., dE/dpc= 0) leads to the equilibrium density 

pc = C^^^ ( A - B M^^^ y^'^, (4.3) 

and stable equilibrium (i.e., d E/dpc > ) exists if and only if the mass is less 
than the Chandrasekhar limit 

M < Mch = ( A / B f'^. (4.4) 

Untill here the conventional theory, involving the postulate of isotropic 
pressure (i.e. x(r) = 0). If we allow for local anisotropy we shall search for pairs 
of functions, p(r) and x(r), fulfilling (2.4) and (3.16), and select amongst them 
the pair which makes the energy a minimum. An exact solution of this problem 
would be quite involved and we shall use a variational (approximate) method 
which we describe in the following. 

We begin by analyzing more closely the hydrostatic equilibrium eq. 
(2.4). Taking (3.16) into account it may be written 

d [p (1-2 x)] /dr = - m po /r^ + 6 p X /r . (4.5) 

A formal integration of (4.5) gives 

X = ( 2 r^ p)-' I' r^ dr [ p' m po / r ' ] • 

(4.6) 



This is an exact equation (within Newtonian gravity) but it is not closed. In 
fact, it is an integro-differential equation because p depends on x (see (3.18)). 
Fixed the function po(r), eq.(4.6) gives the unique solution of (4.5) which is 
regular at the origin ( it begins with x = a r^ + b r'* + ... ). We also need that x is 
regular at the star surface, where the pressure vanishes. An obvious necessary 
condition is that the integral of the right hand side of (4.6) goes to zero when 
r— >R, R being the star radius. This gives, after an integration by parts, 

io^(mpo/r)d^r = 3lo^p d\ ^ 

3 lo^ Pa d\ = lo" ( m po /r ) d\ - 3 lo^ po d\ - 3 Jo'' Pi d\ > . (4.7) 

This relation may be called the virial theorem with local anisotropy, and 
becomes the standard virial theorem when x = 0. The inequality in (4.7) 
follows from the nonegativity of Pa and it is a necessary condition for the 
function x(r), obtained by solving (4.6), being physical (i.e. x g[-1, 1/2] 
everywhere). 

Now we select a family of locally anisotropic star configurations, using 
the change of variables (4.1) and assuming that 0(^) is the solution of the Lane- 

Emden equation of index n=3. Thus for every value of pc we obtain a function 
Po(r) and, hence, we get x(r) by solving eq.(4.6). For any pair {po(r), x(r)} the 
energy of the star becomes 

E =1 { Uo(r) + ui(r) + UA(r,x) - m / r } po d\ . (4.8) 

Performing the integrals we may calculate the energy in terms of the central 
density. The (approximate) equilibrium configuration of the star will be the one 
with minimal energy within the family. All the integrals in (4.8) are 
straightforward except the third one, involving x(r), which cannot be evaluated 
without finding explicitly the solution of eq.(4.6). We get 

E = ( A M - B M^^^) pc M + C M p^"'^^ + J UA(r,x) p(r) d^r . (4.9) 



We may avoid the calculation of the last integral using the virial theorem (4.7). 
In fact, performing the integrals in (4.7), except the one involving pa, we get 

( B M^^^ -AM) pc^^^ + C M pc"^^^ = 3 1 pA(r,x) d^r > 0. 

(4.10) 

On the other hand, a comparison of (3.17) and (3.18) shows that ua p = pa, 
which allows us to get, from (4.9) and (4.10) 

E =-K'M + 2CM pc"^'^ if E>Ei, E = Ei otherwise , (4.11) 

where Ei is given by (4.2). For any mass above the Chandrasekhar limit (4.4), 
the inequality (4.10) (i. e. E > Ei ) is obviously fulfilled (remember that A, B 
and C are positive), but the energy (4.11) has no minimum for finite Pc , it 
always decreases with increasing central density. This means that there are no 
equilibrium configurations with anisotropic pressure and the star is unstable 
against collapse. However, if the mass is below the limit (4.4), the minimum of 
(4.11) compatible with (4.10) corresponds to the density (4.3) and the equality 
sign E = Ei , that is local isotropy. We see that our method agrees with the 
conventional one in this case, and we predict stable equilibrium with Pc given 
by (4.3) below the Chandrasekhar limit (4.4). 

We must point out that both our calculation (leading to (4.11)) and the 
conventional one (leading to (4.2)) rest upon variational approximations, i.e. 
minimizing the energy only within a restricted family of configurations. But 
our accuracy is superior because we select a family of exact solutions of the 
hydrostatic equilibrium eq.(2.4), whilst the conventional treatment selects a 
family of approximate solutions. Consequently, above the Chandrasekhar limit 
our treatment is more correct. But below that limit, where we do not find 
equilibrium configurations with local anisotropy, the conventional treatment 
provides probably a good approximation. Although rigorous results (for some 
results see Ref.5) and/or more exact (numerical) calculations would be 
Wellcome, we may conclude that allowing for local anisotropy in Newtonian 



white dwarfs gives essentially no new result. For any mass below the limit 
(4.4) the equilibrium configuration is locally isotropic, and there are no 
equilibrium configurations with mass above that limit. Newtonian gravity 
predicts the collapse of white dwarfs with mass above the Chandrasekhar limit. 
In the next two sections we show that the conclusion changes dramatically in 
post-Newtonian gravity. 

5. LOCAL ANISOTROPYIN POST-NEWTONIAN APPROXIMATION 

Although the hydrostatic equilibrium equation with local anisotropy is well 
known [1,6], we shall rederive it for the sake of clarity. Solutions of the 
Einstein equations with spherical symmetry may be most easily obtained using 
curvature coordinates. The metric is 

ds^ = e" dr^ + r^ (de^ + sin^e d^^) - dt^. 

(5.1) 

Here we shall consider only static fields, where a and P depend on the single 

2 2 

coordinate r. Using the notation a' for da/dr and a" for d (X/dr , and similar 
for p, Einstein's equations may be written 

-87up, = Gii = r-'[l-e-«(l+rP0], 

- 87U pt = G22= G33 =1/4 e"" [-2 p" - p'V a' p' + 2 r"^ (a'-p')] , 

87tp = G44 = r"^[l -e""( 1 -raO]. (5.2) 

Obtaining the condition of hydrostatic equilibrium from eq.(5.2) is 
straightforward. The third equation allows getting a in terms of the density p. 
Hence the first equation gives P in terms of p and pr. If these expressions for a 
and P are put into the second eq.(5.2) we obtain 



dpr/dr + 2 (pr- pt) /r = - m p/r^ [1 + p/p ] [1+ 47t r^ p/m ] [1- 2 m /r ]"^ 



(5.3) 



m(r) 



47U p dr. 



Obviously this equation becomes the conventional equation of hydrostatic 
equilibrium if we assume local isotropy (pt = pr = p). 

Einstein's eqs.(5.2) do not imply local isotropy, but only a relation 
between anisotropy of pressure and anisotropy of the curvature of space, 
namely 

87C (pt - Pr) = R22 - Rii , R33 = R22 , (5.4) 

where R^v are components of the Ricci tensor, the last equality being a 
consequence of the spherical symmetry. The connection between anisotropy of 
matter pressure and anisotropy of space curvature leads to the conjecture that 
local anisotropy may be more natural and relevant in general relativity than in 
Newtonian gravity. Another argument is that the effective gravitational force 
depends only on position in Newtonian theory, but depends on the velocity 
vector in general relativity (compare the right hand sides of (2.4) and (5.3), 
taking into account the relation (2.3) between stress tensor and linear 
momentum). This conjecture is supported by the calculations reported in the 
present article. 

As is well known, in order to solve the hydrostatic equilibrium eq.(5.3) 
we need, in addition to boundary conditions, a relation between p, pr and pt . 
This should be obtained from the equation of state which, in general relativity, 
is usually written in the form of a relation between the energy (or mass) 
density, p(r), and the number density of baryons. In the study of white dwarfs 
in post-Newtonian (PN) approximation it is more convenient to use the mass 
density po (2.5) rather than the baryon density. The total energy (or mass) 
density will be 



p = po (1 + u) , 



(5.5) 



where u is given by (3.17). 

In PN approximation the hydrostatic equilibrium may be obtained from 
eq.(5.3) by neglecting terms of second order in p/p and m/r. Using (3.16) we 
get 

d[ p (l-2x)]/dr = -mp /r^ - m p (l-2x)/r^- 47t r p (l-2x) - 2 m^p/r^ + 6p x/r. 
(5.6) 

Similarly as we obtained eq.(4.6) from (4.5), a formal integration of (5.6) gives 

x(r) = ( 2 r^ p)"^ e"^^'^ Jo' e^*"''^ r'^dr'[dp/dr'+ mp /r'^+m p/r'^+47tr'p p +2 m^/r'^] 

=(2r'pf Io'"r''dr'{[l+I(r')- I(r)][mp /r'-3p] +2 mVr'']}, I(r)= Ij r'-'d(mrO, 
(5.7) 

where, in the final expression for x(r), we have neglected terms of second 
order, to be consistent with the PN approximation, and performed several 
integrations by parts. The condition that x is finite at the star surface leads, 
after an integration by parts, to the virial theorem in the PN approximation with 
local anisotropy: 

J d^r { [1 - L"^ r' ^ d(mrO] [ mp /r - 3 p ] + 2 m^ pjr^} = . (5.8) 

where we have retained only terms up to first order in p/p and m/r, and 
substituted po for p in the last term as is consistent with the PN approximation. 

The problem of equilibrium closely parallels the Newtonian one and 
may be solved as follows. Taking into account that the radius R of a star is 
defined by the condition p(R) = 0, and therefore pr(R) = Pt(R) = 0, the total 
mass Mtot and the number of baryons N are given by 

Mtot = 47U lo'^ p r^ dr , N = 47U L'^ n [1- 2 m /r]"^^^ r^ dr . (5.9) 



In PN approximation it is convenient to define a mass M that includes only the 
rest mass of nuclei and electrons, in the form 



M= 47cJo' 



[l-2m/r]"^^Vor^dr. 



(5.10) 



The energy of the star is defined as the difference ( Mtot - M) c^. We get (c = 1) 



where, as usual, we have neglected terms of second order. The condition of 
equilibrium with local anisotropy is that the energy E (5.11) is a minimum for 
fixed M (5.10) amongst configurations fulfilling the hydrostatic equilibrium eq. 
(5.6). The practical procedure is described in the next section. 

6. WHITE DWARFS WITH LOCAL ANISOTROPY IN PN THEORY 

We summarize the conventional theory (i.e. postulating local isotropy from the 
very beginning) of white dwarfs in the PN approximation [4]. We should 
evaluate (5.11), using the expression for u given in (3.17) with x=0, perform 
the change of variables (4.1) and assume that 0(^) is the Lane-Emden function 
of index n=3. But there is a difficulty because it is the mass M , given by 
(5.10), rather than the volume integral of po (4.6) which should be fixed in the 
variational approach. A method to solve the problem is to change from the 
radial coordinate r to the new one r' related to the first by 



dV being the invariant volume element. Then it is straightforward, although 
lengthy [4] , to get 



Ei =(AM - BM'") pc" - K'M + C M pc""' - [DiM'" + D2(A/B) M^''] pc 



E = 47U po r^ dr (u - m/r -3/2 mV ) , 



(5.11) 



471 r'^ dr' = dV = 471 [1- 2 m /r]"^^\^ dr , 



(6.1) 



(6.2) 



The values of the constants are given in Ref.4, except Di=0.02543, 
D2=0.89285. The condition that (6.2), as a function of pc , has a minimum 
implies that stable equilibrium is possible only if the mass is below the limit 

M < {PJBf^ [1 - 2 (A C B"'*)^'^] , D = Di+ D2 , 

(6.3) 

which is very close to the Chandrasekhar limit (4.4). 

We pass to the theory with allowance for local anisotropy. We should 

select po(r) as in (4.1) and get x(r) from (5.7). The energy is obtained adding 
to Ei (6.2) the integral of the last term of (2.17), but we may simplify a lot the 
calculation by subtracting (5.8) from (5.11). This gives, after some algebra, 

E =1 d\{poUi -3pi +3Kpo'^^^ m/r -2 pomV +[mpo/r - 3Kpo^^^ ] \S" 
'd(mr')]}.(6.4) 

An additional advantage of (6.4) as compared with (5.11) is that now all terms 
are of first order and the change (6.1) is not necessary. The integrals are 
straightforward using (4.1) and we get 

E = 2 CM pc"^^^+ F M"^Po^'^ if E > Ei , F = Fi (A/B) M"^^^- F2 , (6.5) 

where Fi= 1.2483, F2 = 0.3305. We remark that Fi - F2 = Di + D2 (see 
eq.(6.2)). 

It may be realized that the energy as a function of pc has a minimum for 
any M. However for low M the minimum does not belong to the range of 
central densities pc allowed by the inequality E > Ei , which means that the 
equilibrium configuration possess local isotropy. Consequently, the standard 
theory applies in this range of masses (see comments at the end of section 4). 
For large M there are equilibrium configurations with local anisotropy, which 
we study in more detail in the following. 

The transition from local isotropy to anisotropy happens for a mass 
quite close to the Chandrasekhar limit. (The transition mass corresponds to the 



simultaneous fuUfilement of two conditions, namely that the central density pc 
minimizes E and that E = Ei). For any mass above this limit the equilibrium 
configuration of the star possess anisotropic pressure. The most interesting 

parameters in this case are the central density (obtained by minimization of 
eq.(6.5)) and the radius (obtained from eq.(4.1)). We get 

pc= C M-^'\Fi(Am)-F2M^'Y, R = 2.287 C"^^^ M^^^ [Fi(A/B) -F2 M^^^ ]^^^ (6.6) 

We see that, for stars with the same A (i.e. the same chemical composition [4]) 
the central density decreases with increasing mass above the Chandrasekhar 
limit. The radius also increases but more slowly than the mass, so that the 
gravitational red shift at the surface, M/R, increases with mass. 



7. DISCUSSION 



I propose that equilibrium of a star requires that: 1) the general relativistic 
equation of hydrostatic equilibrium, eq.(5.3), is fulfilled, and 2) the total energy 
is a minimum, for a given number of baryons, with respect to two possible 
variations. The first is the spatial distribution of matter, given in spherical 
symmetry by the function p(r). The second is the distribution of linear 
momentum at every point in the interior of the star, given by x(r) in our 
approach. We remember that both the mass and the linear momentum 
contribute to gravity in the general theory of relativity. In the standard 
treatments the variations of the second kind are not allowed, that is local 
isotropy is imposed as an unjustified postulate. 

I have shown that there are equilibrium configurations of white dwarfs 
above the Chandrasekhar limit. Furthermore the equilibrium corresponds to a 
minimum of the energy with respect to neighbour configurations, which will 
give stable equilibrium if we assume that in any slow oscillation of the star, 
eq.(5.3) holds approximately at all times. The question whether there is also 
stability against rapid oscillations should be investigated more thoroughly. The 
central density predicted for the high mass white dwarfs with local anisotropy 



is smaller than the central density of stars with mass close to the 
Chandrasekhar limit, but it should require a more detailed study the question 
whether that density is low enough to make the stars stable against 
neutronization. If the answer were in the affirmative there might exist stars 
with several solar masses sustained by the pressure of the degenerate electron 
gas. If these stars have become cold enough that they do not radiate, they could 
contribute to the missing mass in galaxies. 

Our study involved essentially two approximations: post-Newtonian 
gravity and small local anisotropy. This has constrained us to deal with white 
dwarfs having masses close to the Chandrasekhar limit. For masses much 
higher a numerical solution of the exact equilibrium equation (5.3) would be 
necessary in order to know whether there are locally anisotropic equilibrium 
configurations. 

The theory of equilibrium with local anisotropy here developped may 
be applied to other systems of astrophysical interest. Supermassive stars with 
local anisotropy are studied elsewhere [7]. A numerical study of neutron stars 
with local anisotropy would be also interesting. Finally I shall comment on the 
implications of the theory for main sequence stars, in particular the Sun. 

The standard model of the Sun involves the "a priori" asssumption of 
isotropic pressure at every point [8]. But it might exist a model of the Sun with 
local anisotropy which is more stable than the standard model and still is 
compatible with all structure equations plus the observational data of mass, 
radius, luminosity, etc. If this is the case, the results of the present paper 
suggest that such a locally anisotropic model would have a smaller central 
density and, therefore, smaller central temperature. This would lead to a 
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reduced production rate of Be and B neutrinos without changing the 
production of p-p neutrinos. In this way, allowing for a local anisotropy in the 
Sun might solve the solar neutrino problem. 
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